Quantum versus classical effects in two-photon speckle patterns 
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We discuss quantum and classical aspects of two-photon interference in light transmission through 
disordered media. We show that disorder is the main factor that suppresses the interference, what- 
ever the quantum state of the incident light. Secondarily, the two-photon interference is affected by 
the quantum nature of light (i.e., the well-defined number of photons in the two-photon entangled 
and Fock states as compared to the coherent state). And finally, entanglement is a resource that 
allows to prepare two-photon states with special symmetries with respect to the interchange of the 
photons and, in particular, the states with bosonic and fermionic symmetries. The two-photon 
interference is more robust for the latter states and its sign can be inverted for the fermionic state. 

PACS numbers: 42.25.Dd, 42.50.Ar, 42.50.Dv 



I. INTRODUCTION 

Multiple scattering of light in random, disordered me- 
dia is a subject of great fundamental and applied impor- 
tance. Most of the studies in this field focus on classi- 
cal phenomena that do not require quantization of the 
electromagnetic field to be understood [TH3]- Light is 
even commonly categorized as a "classical wave" [2J, in 
line with sound and elastic waves, as opposed to "quan- 
tum waves" describing electrons in disordered solids and 
matter waves in cold atomic systems. Meanwhile, even 
if states of light exist that mimic classical behavior very 
closely, light is probably the "most quantum" of all waves 
because it provides an unprecedented freedom in con- 
trolling and measuring its quantum state [3] . Indeed, the 
states of light that do not allow for a classical description 
(e.g., single and entangled photons, squeezed light, etc.) 
can nowadays be generated almost at will [SJ |S] , open- 
ing very promising perspectives for communication and 
information processing applications, including quantum 
computation [7]. 

Multiple scattering of non-classical light in random me- 
dia attracts the attention of physicists since the pioneer- 
ing paper by Beenakker [5] who introduced a convenient 
formalism of input-output relations in this field. In par- 
ticular, propagation of pairs of entangled photons in ran- 
dom media has been studied both theoretically jHHTUj and 
experimentally |111 112) in recent years. Experimentally, 
entanglement can be studied and characterized by coin- 
cidence measurements: one measures the probability P 2 
that two detectors each detect a photon simultaneously 
p~3|. [14] . When two photons pass through some (sim- 
ple) optical system (e.g., a beam-splitter, a delay line, 
or both in sequence), P 2 is affected by quantum interfer- 
ences between different propagation paths leading to the 
same measurement outcome. Entanglement of photon 
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states modifies the result of these interferences and can 
be read from the dependence of P 2 on the parameters of 
the optical system. Replacing the simple optical system 
by a random medium and assuming that the coincidence 
rate P2 is measured as a function of the positions ri and 
r 2 of the two detectors, we arrive at the concept of the 
two-photon speckle pattern P 2 (ri,r 2 ) [9 j. This concept 
can be generalized to non-stationary (e.g., pulsed) light: 
the two-photon speckle pattern becomes time-dependent 
P 2 = ^2(1*1, t\] r 2 , t 2 ). For two independent photons, P 2 
factorizes: P 2 (l,2) = Pi(l)Pi(2), where we abbreviate 
i = {rifti}, i = 1,2, and Pi(i) is the probability to de- 
tect a photon at a position at a time t j . It describes the 
usual, one-photon speckle pattern and it is proportional 
to the intensity of light at r$ . 

It is important to realize that P 2 defined above corre- 
sponds to a single realization of the random medium. It 
is therefore a random quantity and fluctuates from one 
realization of disorder to another. To obtain a determin- 
istic quantity, it is therefore natural to average P 2 over an 
ensemble of realizations of the random medium |15) . This 
ensemble- averaged quantity was studied in Refs. [9HT2]. 
We will denote the ensemble average over random realiza- 
tions of disorder in the medium by a horizontal bar • ■ •. 
It should be distinguished from the quantum-mechanical 
expectation value denoted by (•••). Note that even for 
two independent photons, P 2 (l, 2) = P 1 (1)P 1 (2) does not 
factorize into a product of Pi's (in contrast to the un- 
averaged P 2 ) because P\{i) can have nontrivial (classical) 
correlations in both space and time [DE]. Therefore, the 
ensemble-averaged two-photon speckle P 2 (l,2) combines 
properties due to the quantum nature of the incident light 
and those arising from the classical correlations between 
photons at two different positions (or times). This obser- 
vation is the starting point of the analysis that we develop 
in the present paper. In the attempt to separate quan- 
tum effects from classical ones, we find that properties of 
two-photon speckles are conditioned by four distinct phe- 
nomena. First, the indistinguishability of photons plays 
a very important role. It does not require entanglement 
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and can be present in both quantum non-entangled states 
(e.g., the two-photon Fock state) and the classical (i.e., 
coherent) state. Second, the quantum nature of the inci- 
dent light is important. The quantumness of two-photon 
states, either entangled or not, is expressed by the fact 
that their number of photons (n — 2) is known and con- 
served. In contrast, the coherent state does not have 
a well defined number of photons associated to it. And 
third, the classical correlations between photons, induced 
by the fact that they propagate through the same ran- 
dom medium, turn out to be the most important factor. 
They determine the overall behavior of 2) whatever 
the state of the incident light. As a result, P2 calculated 
for two-photon entangled and non-entangled states, as 
well as for the classical light (coherent state) look very 
similar, with, however, one notable exception. Entangle- 
ment allows us to modify the symmetry of the state with 
respect to the exchange of the two photons [12j [16] . Re- 
sults that are qualitatively different from those for Fock 
and coherent states are obtained for the antisymmetric 
state that changes sign upon this exchange. Photons be- 
have as if they were fermions in this case. This behavior 
cannot be reproduced with non-entangled states and can 
be thus used as a hallmark of entanglement. 



II. THE MODEL 
A. Random medium and input-output relations 

Consider a three-dimensional slab of elastically scatter- 
ing random medium (no absorption or gain), perpendicu- 
lar to the z axis and having thickness L and cross-section 
A ^> L 2 (see Fig. l]) . The incident light is multiply scat- 
tered by the random heterogeneities of the medium be- 
fore reaching one of the two detectors located in the far 
field. In the multiple scattering regime, the mean free 
path I due to disorder is much less than the slab thickness 
L, £ -C L. On the other hand, the disorder is considered 
to be weak, so that £ Ao, with Ao being the central 
wavelength of the incoming light. This corresponds to 
the diffuse regime in multiple scattering pQ. 

Monochromatic light at a frequency uj incident on the 
slab can be decomposed over the basis of plane waves hav- 
ing wave vectors k(w) = {k±,k z }, with kj^ = {k x ,k y }, 
in one of two orthogonal polarization states that we will 
denote by 'o' and 'e'. This notation refers to ordinary 
and extraordinary polarized waves in a bircfrigcnt non- 
linear crystal where, as we will explain below, the two 
waves can be generated (see Fig. [I]) . For a slab of surface 
A, the number of modes in this basis is 2 x N(uj) — 
2 x k(uj) 2 A J '2ir . The same representation is valid for 
the light leaving the slab. From here on we will assume 
that N can be assumed constant (i.e., independent of 
uj) within the frequency bandwidth of the incident light. 
Input-output relations give the photon annihilation oper- 
ators Cbiiw) associated with the outgoing modes i in terms 
of the annihilation operators d a (uj) associated with the 
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FIG. 1. Sketch of the experimental situation. A birefrigent 
nonlinear crystal is pumped by laser pulses of central fre- 
quency 2d) and bandwidth a. Two collinear beams of orthog- 
onal polarizations 'o' (ordinary) and 'e' (extraordinary) and 
frequencies uj\ and u)2, such that Wi + 012 = 2w, are obtained 
as a result of spontaneous parametric down-conversion in the 
crystal. These beams are incident on a slab of disordered 
medium, with the 'e' beam delayed by a time r. Light is 
multiply scattered inside the slab and the transmitted light 
in modes i and j is detected by two photodetectors in the 
far field. The multiple scattering regime in ensured by the 
requirement <<L, with £ the mean free path due to disorder 
and L the thickness of the slab. 



incoming modes a |S]: 

N 2N 

di (uj) = ^ *o» M a a (uj) + ^ r 0i M °>P M • C 1 ) 

a=l 0=N+1 

For the outgoing mode i on the right of the slab (as in 
Fig.[lj), the first sum runs over the incoming modes a inci- 
dent from the left and transmitted through the slab with 
transmission coefficients t a i(ui), whereas the second sum 
runs over the incoming modes j3 incident from the right 
and reflected with reflection coefficients rpi(uj). The co- 
efficients t a i(u>) and rpi{uj) form the (unitary) scattering 
matrix S of the slab. Operators di and a a obey the usual 
bosonic commutation relations: [di(u>i), aj(u>2)] = and 
[<Zj(a>i), dj(uj2) ] = 8ij5(uii — oj 2 ), and the same for a a 
and dp. 

As far as L ^> £, N 3> 1, and we are not interested 
in quantities that involve summations over all input or 
output modes (like, e.g., the conductance or the total 
transmission), the unitarity constraint on the matrix S 
can be relaxed. Then, in the diffuse regime of scattering, 
the transmission coefficients t a i(uS) can be assumed to 
be independent identically distributed random variables 
with circular Gaussian statistics, zero mean, and some 
frequency correlation function C(Au>): t a i(uS)t a i j{uj')* — 
TS aa 'SijC(uj — uj'). Here T = |i Qi (a;)| 2 is the frequency- 
independent, average intensity transmission coefficient. 



B. Quantum states of light 

The focus of this paper is on the two-photon entangled 
state that can be experimentally obtained by the spon- 
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taneous parametric down-conversion (SPDC) in a non- 
linear optical crystal [H] HZ]- To be specific, we will 
restrict our consideration to the collinear type-II SPDC 
process in which two entangled photons have orthogonal 
polarizations ('o' and 'e') and propagate collinearly with 
the pump beam |18j . At low intensity of the pump, the 
two-photon state generated in the SPDC process can be 
written as jTSj 



duiiduj2B (wi,w 2 ) Q-o ( w i)^ (w 2 ) |0), (2) 



where B {u)\, UJ2) — Ka (u>i + ljz) $ (u?i, U2) |19j - 

a (wi + uj 2 ) = exp[-(cji + w 2 - 2iJ) 2 /2a 2 } (3) 

is the pump envelope and 

<E> (wi , 0J2) — sine [v a (cji — uj) + v c (w 2 — uf)] (4) 

follows from the phase matching condition [TB]. Here 
2uJ is the frequency and a is the spectral width of the 
pump beam; Vj = L NL (dkj/du\ u = G - dk p / 'du}\ u = 2 u) /2, 
with j = o, e, quantify the phase mismatch between the 
down-converted photons 'o' and 'e' and the pump 'p' in 
the nonlinear crystal of length Lnl- kj(oj) denotes the 
dispersion relation of the crystal for ordinary (J = o) 
and extraordinary (j = e) polarized light, as well as for 
the pump (j = p). The constant K follows from the 
normalization condition (V'entlV'ent) — \K\ 2 Tr 3 / 2 a/\r]_ \ — 
1. 

The entangled nature of the state ^ is manifest in 
the fact that the function B(u>i,uj2) is not factorizable 
in a product of two functions of w\ and w 2 , respectively. 
The two photons in the state (2j are thus entangled in 
frequency. A state similar to (12) was recently considered 
in Ref. [TU]. In contrast to that work, however, our Eq. 
([2]) describes a state that is not symmetric with respect 
to the exchange of frequencies of ordinary and extraor- 
dinary photons: B(u>i,W2) ^ -B(w 2 ,u;i). This property 
stems from the different dispersion relations for the two 
polarization states in the nonlinear crystal and can, as we 
will see, have important consequences for the probability 
P 2 of simultaneous photon detection by two photodctec- 
tors because the two photons become more and more 
distinguishable as the bandwidth of the pump pulse a is 
increased. We illustrate this in Fig. [2] where we show the 
spectral width Aw, defined as the full width of the spec- 
trum at half-maximum height, of ordinary and extraordi- 
nary beams, normalized by their value Aw cw = 2.78/|?y_| 
in the monochromatic limit a —> 0. To be able to esti- 
mate this quantity analytically, we replace the square of 
the sine function in Eq. Q by a Gaussian exp(— x 2 /2.79) 
with the same width at half-maximum. We obtain 



2|i/„ 



\V- 



\ 



ln(2) 



2.79 



K, e Aw cw ) 2 V A ^c 



,(5) 



where ?7_ = v Q — v c measures the typical delay between 
the two entangled photons when they leave the nonlin- 
ear crystal where they are generated. If the pump is 
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FIG. 2. Spectral widths of the ordinary and extraordinary 
beams of photons in the state described by Eq. |2|, as func- 
tions of the bandwidth of the pump beam a, for typical pa- 
rameters of a BBO nonlinear crystal [20]. All quantities are 
normalized by the spectral width Aa> cw of the two beams in 
the monochromatic limit a = 0. 



monochromatic, \r\-\ can be considered as the coherence 
time of the 'biphoton' [17] . We see from Eq. ^ and Fig. 
[2] that the spectral widths of the two beams become dif- 
ferent when a increases, making distinguishable the pho- 
tons originating from different beams after transmission 
through a disordered medium. 

In order to understand the role of entanglement in the 
final results for P 2 , we will compare them to calcula- 
tions performed for non-entangled states of light. The 
first non-entangled state that we will consider is the two- 
photon Fock state |"0Fock) given by Eq. ^ with fac- 
torizable B(uji,uj2) — B{uo 1 )B{uj2). We will consider 
B(ui) = exp[-(ui - w) 2 /2A 2 ]/(V^A) 1 /2 ; w here A is the 
bandwidth. This state has the same number of photons 
n = 2 as the entangled state |2]) but does not imply any 
quantum correlations between them. However, the two 
photons have the same spectrum, so that they are in- 
distinguishable after transmission through a disordered 
medium. 

Finally, the simplest and the most widespread non- 
entangled state is the two-mode coherent state 



|^coh)H{a(u;)}) |{a(u;)}) e 



(6) 



with a (w) = exp[-(a;-w) 2 /2A 2 ]/(0i : A) 1 / 2 [5T]. It is an 
eigenstate of the photon annihilation operators a (uj) and 
a c (uj): a ,o(w)|V'coh) = a ( w ) IVWi)- This state describes 
two pulses with orthogonal polarizations. The funda- 
mental difference between the coherent state Q and the 
previously introduced two-photon states is that the for- 
mer is not an eigenstate of the photon number operator 
h. The number of photons is therefore not a good quan- 
tum number in this state. However, similarly to the Fock 
state, the photons in the state ^ have the same spectral 
properties and hence cannot be distinguished after being 
transmitted through a disordered medium. 
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III. AVERAGE PHOTOCOUNT RATE 

The entangled, Fock and coherent states introduced 
above have (on average for the coherent state) one photon 
in each polarization state. This can be confirmed by 
calculating the quantum expectation value (h 0;B ) of the 
photon number operator |21j 



r AT/2 

/I,,., = I dta\ c (t)a . c (t), 

-AT/2 



where hi(t) is the inverse Fourier transform of di(ui). This 
quantity can be measured in an experiment by counting 
photons arriving at a polarization-discriminating pho- 
todetector placed in front of the photon source during 
a sufficiently long time AT. 

To start with and to illustrate our general calculation 
scheme, let us compute the average number of photons 
transmitted into a mode i when the random medium is 
illuminated by light in one of the three states |'0ont), 
I^Fock) or IV'coh) defined in the previous section. The 
operator corresponding to this quantity is |21) 



rAT/2 

a, : ; dt a\(t)ai(t). 

'-AT/2 



(8) 



To perform the calculation, we express the operator hi (t) 
through the operators a a (t) and ap(t) using the input- 
output relations ([I]), use commutation relations to com- 
pute the quantum expectation values and the Gaussian 
statistics of t a i to average over realizations of disorder. 
The result is the same for the three quantum states that 
we study: 



(hi) = ((h ) + (h c ))T = 2T. 



(9) 



We thus observe that the average photon number is not 
sensitive to the quantum state of light, — a well-known 
and quite general fact in quantum optics. 

Given the normalization of the input states and the 
fact that T <C 1 in the diffuse regime, (hi) is equal to the 
ensemble-averaged probability Pi that a photodetector 
measures a photon in the outgoing mode i. It also gives 
the photon counting rate (in units of photons per pulse) if 
a sequence of independent and separated in time pulses in 
quantum states \ip cn t), IV'Fock), or |V>coh) is sent through 
the random medium. 



IV. COINCIDENCE RATE AND PHOTON 
CORRELATIONS 



Two different but related quantities can be used to 
characterize correlations between numbers of photons 
transmitted into two modes i and j: the correlation func- 
tion Cij — (hihj) and the probability -P2 («,:/) to detect 
a photon in each of the modes i and j. The later also 
corresponds to the rate of photon coincidence counts if a 
sequence of pulses in the same quantum state is sent into 



the medium. The two quantities are related through a 
relation that we derive in Appendix [A] 



P2(i,j) 



Cij — 5ij (hi) 



1 



1 



(10) 



where : • • • : denote the normal ordering of operators. 
It is worthwhile to note that P 2 (i,j) was identified with 
(: hihj :) in some of the recent literature |10j without 



(7) distinguishing the cases i ^ j and i = j 



As we see from Eq. (10 1, both P 2 (i,j) and CV, can 
be found from the normally ordered correlation function 
(: hihj :). The average of the latter over disorder can be 
computed for the three quantum states that we defined 
in Sec. II B (see Appendix [Tj] for details): 



2T 



1 



x B{u 2 ,u 1 )\C(u 1 -u 2 )\ 2 e- i ^- u '^ 

niftj : >Fock = 2 ^ 1 + *• 

x |BM| 2 |C( Wl -u 2 )\ 2 e- 1 ^^ ^ 
hihA :)__. = 4T 2 1 + 8, 



du)\dui 2 B* (lui , 012) 

(11) 

duJiduj2\B(uii)\ 2 

(12) 

du>idu>2\a(uji)\ 2 



x |a(w 2 )| 2 |C(w 1 - w 2 )| 2 cos 2 ([ Wl - wa]r/2) 



(13) 



In the multiple scattering regime, the quantities given 
by Eqs. (|r"Lj)-(fl3|) will be very small and one can doubt 
their interest for real experiments. Indeed, typically we 
have T - t/NL with L/t > 10. The number of trans- 
verse modes for a slab of area A = 1 mm 2 will be 
= fc 2 yl/27r ~ 10 7 at optical frequencies, leading to 

a prefactor T ~ 10~ 16 in Eqs. Neverthe- 
less, photon coincidence measurements in transmission 
through a disordered medium were realized for coherent 
and quasi-chaotic states [22j. The measurement can be 
optimized by noting that Eqs. |TT])-( 13 ) do not depend 
on the particular choice of modes i and j but are only 
sensitive to the fact that i = j or i ^ j. One can thus 
count photons in an arbitrary large number M > 2 of 
modes and then average over the results obtained for all 
pairs (i,j) of modes. Each mode should, however, be 
addressed individually, which will require M photode- 
tectors. One can also think about more sophisticated ex- 
perimental setups to detect photon coincidences without 
knowing the precise mode to which the photons belong, 
like the one exploiting the two-photon absorption [23] - 
From the theoretical point of view, we can avoid working 
with too small quantities by normalizing Eqs. (11 )— ( 13 ) 



by T . We thus define the normalized photocount coin- 
cidence rate R as 



Ri 



1 



1 



(14) 
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An expression similar in structure to Eq. (11) was de- 



rived for the coincidence rate in Ref. [10] . though for a 
slightly different entangled state. In addition, the au- 
thors of Ref. [TO] have taken into account the lowest- 
order corrections to the Gaussian model for transmis- 
sion coefficients t a i adopted here, and found additional 
contributions to the coincidence rate. These contribu- 
tions, however, turned out to be of order l/g, with 
g = N£/L ^> 1, and would be negligible under condi- 
tions of diffuse scattering that we consider here. In the 
absence of these terms, we find from Eqs. ( 11 )-( 14 1 that 
for i ^ j, R cnt — i?Fock = 2 and i? C oh = 4 independent 
of any parameters. The latter result corresponds to the 
total absence of any correlation between photon numbers 

— 2 

in two different modes i =^ j: (: hifij :) = 4T = (hi)(hj), 
whereas the former indicates that there are negative cor- 

2 

relations: (: hihj :) = 2T < (hi)(hj). Negative cor- 
relations in transmission of non-classical light through 
disordered media were already predicted theoretically in 
Refs. 24-26 and observed experimentally in Ref. [S7J for 
squeezed light. 

For i ^ j, the parameter-dependent corrections to Eqs. 
(fTTb— (fl3]) due to the non-Gaussian statistics of t c 



and 

correlations between them would be of order 1/j C 1 
[TU] , This is much smaller than the variations of the co- 
incidence rates by 50% or so for i = j. The latter should 
be therefore much easier accessible experimentally. For 
this reason, in the rest of this paper we will study the 
case i = j and compute Ra which describes the prob- 
ability for two photons to be found in the same out- 
going mode. We will omit the subscript Hi' of Ra to 
lighten the notation. At this point, the lacking ingredi- 
ent in Eqs. ( fll] ), (12) and (13) is the frequency correla- 
tion function C(AtH). We will use two different models 
for it. A simple model I, C(Ato) — exp(— | Aw|/f2 corr ), 
with f2 corr the correlation frequency, will facilitate in- 



tegrations in Eqs. ( 11 )— ( 13 1 and thus will allow for a 
number of analytic results. A more realistic model II, 
C(Aui) = v /-iAw/fi T h/sinh( v /-iAw/fi T h) [28], with 
fixh the Thouless frequency, will be used to obtain pre- 
dictions that can be experimentally verified in a diffusely 
scattering disordered medium. We note, however, that 
both correlation function given above decay fast with Au; 
and hence the results will be similar for both models. 



A. Entangled state 

For the entangled state |2]) and the exponential model 
I for the correlation function C(Auj) of transmission co- 
efficients, we obtain from Eq. p"T] ): 

R ent (t,s,w) 

If 1 rs 1 

= 1 + — 7= / dxf (t,w,x Erf -(1-M) ,(15) 

where f(t, w, x) — 2w/[A + w 2 (x + t) 2 ] and the dimen- 
sionless variables are t = t/t}—, s = ari+ and w — 
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Time delay r/?7_ 



FIG. 3. Photocount coincidence rate for two detectors in 
the same outgoing mode and the two-photon entangled state 
([2| at the input of a disordered medium with frequency- 
independent transmission coefficients t a i, as a function of de- 
lay t between the two photons. Three curves correspond to 
different values of the parameter s — ar/+ measuring the spec- 
tral width of the pulse that pumps the nonlinear crystal where 
the photons are generated. An infinitely large value of Thou- 
less frequency is assumed. 



I^corr^-I, with r)± = v ± v e . 

For Qcorr ^ 1/1*7-1) ^cm( w ) can be assumed frequency- 
independent within the bandwidth l/|?y_| around u and 
Eq. ( 15 ) reduces to 



Rent(t, S, OO) 



1 



Erf [§ (l-|*l)] 




1*1 <1 
1*1 >1 



(16) 



Exactly the same result is obtained for the model II in 
the limit of | S^Th 1 ?— I ~~ ^ oo. Up to a sign in front of the 
second term in the r.h.s., Eq. (16) is equivalent to the 



coincidence rate in the Hong-Ou-Mandel interferometer 
[18] . Instead of a dip in the coincidence rate as a function 
of delay t, we now find a peak, which is explained by the 
fact that we compute coincidences in the same outgoing 
mode and not in two different modes. We show Eq. ( 16 ) 
in Fig. [3] for 3 values of s, corresponding to 3 different 
bandwidths a of the pump pulse at fixed rj + . Increasing 
a suppresses the two-photon interference effect in R cnt . 
This is due to the loss of indistinguishability of the two 
photons, which, in its turn, results from the fact that 
the two entangled photons of the state ([2| do not have 
the same spectrum (see Fig. [2]). As a consequence, the 
two-photon interference between them washes out with 
increase of a. 

The minimum and the maximum values R = 1 and 
2 of the normalized coincidence rate in Fig. [3] can be 
understood without appealing to the entangled nature 
of the quantum state and turn out to be universal for 
all two-photon states. For |r| > \r)-\, the photons are 
well separated in time and are transmitted through the 
disordered medium independently. The probability for 
each photon to end up in a given outgoing mode i is equal 
to T. The probability that both photons are measured 



() 



^2.0 

a 

*1.8 
cd 
+= 

*-< 1.6 

V 

g 1.4 

• i— I 

p| 1.2 



o 
U 



1.0 



0-77+ = 




— r^ThT?- = 00 






0.3 \ 






o.i — 7 






- 2 - 1 


_i i i i 




1 2 



Time delay r /rj- 



FIG. 4. Photocount coincidence rate for two detectors in the 
same outgoing mode and the two-photon entangled state |2| 
at the input of a disordered medium. The pump beam that 
generates the pair of photons in a nonlinear crystal is assumed 
to be monochromatic (a = 0); different values of |fiThT?- 1 cor- 
respond to different disorder strengths and/or sample thick- 
nesses. 



in the same mode i is equal to the product of single- 

2 

photon probabilities: P2(i,j) — T , leading to R = 1. In 
contrast, at t = the two photons cannot be considered 
independent and the two-photon interference leads to a 
larger probability for them to be detected in the same 
mode, thus the peak in Fig. [3] In the ideal case of ar/ + = 
this interference is perfectly constructive and R doubles 
with respect to the case of independent photons. We will 
show in Sec. [V] that destructive two-photon interference 
can lead to R < 1 and even to R = for certain specially 
prepared states. 

Let us now consider the case of arbitrary relation be- 
tween the Thouless frequency Oxh and r]_ . In Fig. |4j we 
plot the coincidence rate for 3 different values of |f2xh?7- 1- 
We observe that the increase of disorder (quantified by 
the decrease of the parameter |f2xhf7-|) has two distinct 
consequences. First, we observe lowering of the maxi- 
mum coincidence rate in the absence of delay (r = 0). 
The disorder thus washes out the two-photon interfer- 
ence. However, at the same time we see that the curve 
-Rental") broadens and the photon coincidence rate grows 



up at large |t| 



> 



n- 



This is a consequence of the fluc- 



tuating time of flight of photons through the medium. 
Typically, the time of flight of a photon through a dis- 
ordered medium fluctuates from very small values (bal- 
listic propagation) to values exceeding (diffusion). 
The medium can, therefore, partially compensate for the 
initial delay r between the two photons and make them 
interfere even for |r| > |f7-|. 



B. Separable state 

As follows from Figs. [3] and |1J the photocount coin- 
cidence rate for two photodetectors behind a disordered 
medium is sensitive to the parameters a and r\- that 
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FIG. 5. Same as Fig.|4]but for the separable, two-photon Fock 
state incident on the disordered medium. The delay time r 
and the Thouless frequency flTh are in units of bandwidth A 
of incident photon beams. 



characterize the state of the two entangled photons inci- 
dent on the medium. But let us check what really changes 
if now we consider the two-photon Fock state IV'Fock) de- 
fined in Sec. II B This state is separable (not entangled), 
but the two photons are still indistinguishable. There- 
fore, the comparison of results obtained for the states 
l^ent) and IV'Fock) should allow revealing the role of en- 
tanglement in the experiment depicted in Fig. [l] For the 
exponentially decaying correlation function C(Aw), we 
obtain from Eqs. (flOl and (fl2|): 



R F ock(t,w) = 1 + e ^ |t|2 ( ™ )2 ]|cos(^ 



Re e" 21 y Erf 



2]t| 
w 



1/2 
V2 



.(17) 



where t = tA and w — J7 corr /A. 

In Fig. [5] we show the coincidence rate i?Fock for the 
more realistic model II. i?Fock is shown as a function of 
normalized time delay at three different values of the nor- 
malized Thouless frequency. It is bounded by 1 from be- 
low and 2 from above for the same reasons as i? e nt (see 
the discussion in Sec. IV A). Figure [5] is to be compared 
with Fig. [4] At weak disorder (large f^xh) the results for 
the entangled and separable two-photon states are quite 
different: the triangular shape of i? ent (r) is replaced by 
a rounded, Gaussian-like profile of JiFockM- In contrast, 
at small J7xh the two states lead to similar results. We 
thus can obtain very similar dependencies R(t) for two- 
photon entangled and separable states, provided that we 
adjust the bandwidth of the latter to values A ~ l/|r?_|. 
Given this observation, we are led to conclude that, at 
least in the presence of sufficiently strong disorder (i.e., 
for |f2xh?7-|j ^Th/A < 1), the impact of entanglement 
on the coincidence rate is quite limited. In an experi- 
ment measuring R(t) behind a random medium, it would 
be difficult (if possible at all) to conclude if the incident 
state was the entangled state IV'cnt) or the separable state 
IV'Fock) when no additional information (e.g., the band- 
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FIG. 6. Same as Figs. [4] and [5] but for the coherent state Q 
incident on the disordered medium. The delay time r and the 
Thouless frequency S7tii are in units of bandwidth A. 



width of the incident photon pulses) is available. 



C. Coherent state 

Because in transmission through a disordered medium 
the two-photon entangled and separable states lead to 
very similar dependencies of the coincidence rate on the 
delay time between the two photons, we would like to 
clarify the common property of these states responsible 
for this result. Clearly, both states \ip on t) and IV'Fock) 
contain exactly two photons and these two photons are 
indistinguishable. But what is more important? To an- 
swer this question, we now consider the coherent state 
iV'coh) in which the number of photons is not well-defined, 
but the photons are still indistinguishable. Assuming the 
model I for the correlation function C(Aw), we can again 
compute the coincidence rate analytically from Eqs. ( 13 ) 
and (Pll: 



R coh (t,w 




, (18) 



where t = tA and w = f2 corr /A. 

In Fig. [6] we show i? co h as a function of delay time 
t for the realistic model II of the correlation function 
C(Auj). Despite the fact that the coherent state exhibits 
neither entanglement, nor a well-defined photon number, 
the overall shape of the coincidence curve is still the same: 
it is a bell-shaped curve with the overall amplitude sup- 
pressed by disorder. However, i? co h reaches larger abso- 
lute values, up to i? co h = 4, in contrast to i? on t and i?Fock 
that are bounded by 2 from above. 

The results for the coherent state can be rederived in 
a completely classical framework, by replacing the op- 



erators a and with complex numbers a and a* and 
ignoring the commutation relations. This provides a 
way to understand some of the aspects of Fig. [6] using 
the well-known facts from the theory of classical wave 
scattering [HI2]- Indeed, we obtain from Eq. ( [Ft) that 
R co h = 2I^ t /(Iat) 2 , where I at is the intensity of the 
transmitted wave, integrated over a time interval AT. 
When t = 0, the transmitted wave a, is a superposition 
of independent waves a\ and resulting from the 
transmission of incident waves with 'o' and 'e' polariza- 
tions, respectively. For weak disorder, both and a^f 1 
(and hence a,) have zero- mean circular Gaussian statis- 
tics. Because the intensity / = |ai| 2 , I at /{Iat) 2 = 2 is 
obtained for the monochromatic illumination (A — > 0). 
This corresponds to the so-called large intensity fluctua- 
tions, well-known for classical speckle patterns [55]. The 
resulting maximum value of i? C oh is 4, as can also be seen 
from Fig. [6] In contrast, when |r| is large, the detections 
of signals due to the incident 'o' and 'e' waves are sepa- 
rated in time. The measured time-integrated intensity is 



a sum of two independent terms: I at — ^at + ^at- This 
yields /^/(Iat) 2 = 3/2 and i? coh = 3 for A -> 0. When 
the spectral width A of the incident waves increases, the 
measured signals suffer from a partial averaging even for 
a single realization of disorder. In the limit of A — > oo, 
the intensity of transmitted light does not fluctuate and 
-Rcoh = 2. 



r(e) 



V. SYMMETRIC AND ANTISYMMETRIC 
ENTANGLED STATES 

As we have shown in the previous section, the entan- 
gled state ^ with the function B(ui , oj 2 ) defined by Eqs. 
([3]) and Q leads to the coincidence rate R measured be- 
hind a random medium that is similar to that for a sep- 
arable state (e.g., compare Figs. [4] and [5]). In an experi- 
ment, if one simply measures R as a function of delay r, 
it might be difficult or even impossible to conclude on the 
presence of entanglement between the two incident pho- 
tons without knowing their other characteristics, such as 
the bandwidth. This might give an impression that en- 
tanglement adds no interesting physics to the optics of 
random media. This impression is, however, completely 
wrong. In order to demonstrate this, let us now consider 
a little bit more sophisticated entangled state, namely, 
the state given once again by Eq. ([2| where we replace 
B(ui,u} 2 ) by 



B fl (wi,wa) =JC[B(cj 1 ,L) 2 ) + e ie B(L> 2 ,u> 1 )}, 
with K, being a normalization constant, 

1 1 

2 X 1 + cos(6»)Erf (a? ]+ /2)^/a7 1+ 



\ic\ 2 



(19) 



(20) 



By properly adjusting the value of 9, we can make this 
state symmetric (9 — 0) or antisymmetric (9 = n) with 
respect to the exchange of the two photons. In other 



words, the state preserves its form with the same (sym- 
metric) or opposite (antisymmetric) sign upon the ex- 
change u>i 4-> uj 2 and can thus also be said to have bosonic 
or fcrmionic symmetry, respectively. An intermediate sit- 
uation with < 8 < tt corresponds to asymmetric states; 
8 = 7r/2, for example, yields a state that leads to the 
same result for the average two-photon coincidence rate 
as the state ^ that we studied in Sec. IV A 

The states corresponding to Eq. ( 19 ) can be prepared 
experimentally 16] and for 8 = 0, tt have the advantage 
of conserving the perfect indistinguishability of the two 
photons whatever the bandwidth of the pump a. This is 
in contrast to the entangled state considered in the pre- 
vious sections, in which the two photons become distin- 
guishable as a is increased because their spectra become 
different (see Fig. |2j. As a result, Eq. (19) leads to a 



much weaker dependence of the coincidence rate on a. 
In addition, the sign of the two-photon interference term 
can be inversed for the antisymmetric state (9 = tt) , thus 
turning the constructive interference into the destructive 
one. 

Because Bg (u>i , LJ2) is a linear combination of B (oji , LJ2) 
and B(uj2,loi), the calculation of the coincidence rate for 
the former can be reduced to the analysis that was per- 
formed in Sec. |IV A For the exponentially decaying cor- 
relation function C(Acj) we obtain 

Rg(t,s,w) = 1 + 2%|/C| 2 J dxf(t,w,x) 

x [I(s > x) + cos(6)J(s,x)] , (21) 



where f(t,w,x) was defined in Sec. IV A The functions 
I{s,x) and J(s,x) are 



I(s,x) 



1 



=Erf 



(i-N) 



J(s, x) = — (1 — |x|) exp 

TT 



(22) 
(23) 



The striking difference between symmetric (8 = 0) and 
antisymmetric (8 = tt) states is demonstrated in Fig. [7] 
We see that for a monochromatic pump (er = 0) and weak 
disorder (|f2Th?/-| — > 00), the symmetric state produces 
the same dependence of the coincidence rate on the delay 
time r as the asymmetric state considered in Sec. |IV A"} In 
contrast, the antisymmetric state (8 = tt) leads to a dip 
in the coincidence rate at r = 0, instead of a peak. This 
is indeed reminiscent of the behavior of a pair of fermions 
that obey the Pauli principle and therefore avoid being 
in the same quantum state. Instead of an increase of 
coincidence rate observed for the state with the bosonic 
symmetry, the fermionic symmetry results in a complete 
suppression of coincidences, making it impossible to find 
two photons in the same outgoing mode. This effect is 
suppressed when the strength of disorder is increased (see 
the pair of lines corresponding to |fiTh*?-| = 0.3 in Fig. 
[7]), but the dip at t = is still clearly visible. 

The role of symmetry in the structure of the entangled 
state becomes obvious when we look at the evolution of 
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FIG. 7. Photocount coincidence rate for two detectors in the 
same outgoing mode and the symmetric (9 = 0, blue lines) 
and antisymmetric (9 = tt, orange lines) entangled states at 
the input of a disordered medium. The pump is monochro- 
matic and the disorder is weak QQt^t]- | — > 00) for solid lines. 
|^Thf?-| = 0.3 for dashed lines. 
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FIG. 8. Photocount coincidence rate in the absence of time 
delay (r = 0) and for |fixh??-| — > 00 (solid lines), 1 (dashed 
lines) and 0.3 (dotted lines), as a function of the normalized 
bandwidth of the pump [<T77+|. 



Re(r) with the bandwidth a of the pump. According to 
Fig. [2j increasing a makes the spectra of the two pho- 
tons different and thus makes the photons distinguish- 
able. This suppresses two-photon interference effects and 
reduces the peak in i? ent (r) (see Fig.[3|. However, the sit- 
uation is quite different for the symmetrized states that 
we study in the present section: in the states correspond- 
ing to 8 = or tt, the two photons remain indistinguish- 
able independent of a. The dependence of Re(r = 0) on 
|<7?/ + | shown in Fig. [8] illustrates this quite convincingly. 
At weak disorder (|rixh?7-| — > 00), Rg(r = 0) does not 
vary with a, whereas at stronger disorder, increasing the 
bandwidth eventually suppresses interferences, though at 
considerably larger scales as compared to the asymmetric 
8 = tt/2 state. 

To illustrate the positive impact of the symmetry of the 
state on the two-photon interference, in Fig. [9] we show 
the dependence of the photocount coincidence rate on the 
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FIG. 9. Comparison of photocount rates for states corre- 
sponding to 9 = (symmetric state, solid lines) and n/2 
(dashed lines), for a particular value of pump bandwidth 
Wv+ \ = 4- 
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FIG. 10. Comparison of photocount coincidence rates cor- 
responding to 6 — 0, jo-77+l = 4 (solid lines) and 8 — n/2, 
1 077+ 1 = (dashed lines). 



delay time for the symmetric state {8 = 0) and the asym- 
metric state corresponding to 8 = n/2, at three different 
disorder strengths. For weak disorder QQ^hV— I — ^ °°)> 
the difference is very important, although at stronger dis- 
order the two results start to be closer, showing that dis- 
order tends to reduce the role of symmetry. To quantify 
the impact of symmetry even further, in Fig. [10] we re- 
plot the solid lines corresponding to 8 — and \<yq + \ = 4 
and compare them with the result corresponding to the 
8 = n/2 state at \crrj + \ = (monochromatic pump). The 
curves of each of the 3 pairs, corresponding to different 
disorder strengths, are very close to each other. This il- 
lustrates that the symmetry of the state can compensate 
for the loss of two-photon interference due to the large 
bandwidth a of the pump pulse. 



VI. SIGNS OF NONCLASSICAL LIGHT 

Let us now clarify which properties of the coincidence 
rate R measured in transmission through a disordered 
medium result from the classical or nonclassical nature 
of incident light. We adopt the following operational def- 
inition [30l[3T]: we call light nonclassical if its photocount 
statistics p(n) cannot be obtained from the semiclassical 
Mandel's formula [P4] 



p(n) 



-exp(-I)V(I)dI 



(24) 



with V(I) > — the probability density of the classical 
variable (intensity) /. Obviously, the probability to de- 
tect 2 photons in a given outgoing mode i readily follows: 
p(2) = 7 2 exp(-/)/2 ~ P/2. As before, the vertical line 
denotes averaging over realizations of disorder that in 
the present context is equivalent to averaging over the 
distribution V(I) and we made use of the fact that in 
transmission through a thick disordered medium, V(I) is 
appreciable only for I <C 1. We see that Eq. (24) im- 



plies that the normalized photocount coincidence rate as 



defined by Eq. (14 1, is given by 
P(2) 



R = = 2 1 + 5P/I 
T 



(25) 



where 51 = I — I is the fluctuation of intensity and we 
used the fact that I = 2T for the states considered in 
this work [see Eq. d9j]. 

In fact, we already used Eq. (25) to discuss the re- 
sults obtained for the coherent state in Sec. IIV CI with 
I — I at- R = 2 was obtained in the absence of intensity 
fluctuations (5P = 0) and R = 4 — for SP/(I) 2 = 1. 
We thus immediately conclude that the results for the 
coherent state can be described by Eq. (25) following 
from the Mandel's formula (24), confirming that, not 
surprisingly, this state remains classical after transmis- 
sion through the medium. In contrast, we saw in Sees. 
|I V A[ [TVBl and |V| that two- photon entangled and separa- 
ble states lead to R < 2. This would require SI 2 /(I) 2 < 
in Eq. ( |25[ ) , which is impossible for any probability distri- 
bution V\I) . Therefore, the results for two-photon states 
cannot be described by Eq. (24). This shows that these 



states remain nonclassical upon transmission through a 
disordered medium. 

Measurement of the absolute value of R in an experi- 
ment may be complicated because it requires proper nor- 
malization of the photocount number which is the raw 
output of the measuring device. The difference between 
classical and nonclassical states of light also shows up in 
the contrast (or visibility) of the coincidence curve R(t) : 



V = 



\R(0) - R(oo)\ 
R(0) + R(oo) ' 



(26) 



This quantity is not sensitive to the normalization of R 
and may be easier to access experimentally. The max- 
imum contrast is reached in the limit of J7tii 00 ■ 



10 



For the asymmetric entangled state pi), the symmet- 
ric entangled state [Eq. (19 1 with 9 = 0], and the Fock 



two-photon state, the maximum contrast is max(14 nt ) = 
max(Ve = o) — niax(VF OC k) = 1/3, whereas for the coher- 
ent state we find max(V r co h) = 1/7. The largest contrast 
max(Vg =w ) = 1 is reached for the antisymmetric entan- 
gled state [Eq. (19 1 with 9 — %]. Therefore, the contrast 



of the two-photon speckle pattern is more than a factor of 
2 larger for the nonclassical light considered in this paper 
than for the classical light, represented by the coherent 
state. 



VII. CONCLUSION 



The purpose of this paper was to discuss quantum and 
classical aspects of two-photon interference behind a dis- 
ordered medium. We consider a particular realistic ex- 
ample of a light pulse of finite bandwidth sent through 
a disordered medium and containing 2 orthogonally po- 
larized photons, of which one can be delayed by an ar- 
bitrary time t. By comparing the rates R(t) of coin- 
cident photon counting in the same outgoing mode in 
transmission through the disordered medium for differ- 
ent quantum states of the incident light (two-photon en- 
tangled and Fock states, coherent state), we convincingly 
demonstrate that disorder is the main factor that shapes 
the curve R(t). Moreover, if no additional information 
about the incident pulse is available, it is barely possi- 
ble to see any sign of entanglement in R(t) for a typi- 
cal two-photon entangled state generated by the collinear 
type-II SPDC. Almost identical R(t) curves can be ob- 
tained for an entangled and separable (Fock) two-photon 
states, provided that the bandwidth of the latter is ad- 
justed. However, the result obtained for the coherent 
state is different, which highlights the second important 
factor that plays a role in the two-photon interference: 
the quantum nature of light. Here it is manifest in the 
fact that the number of photons is well-defined in the 
two-photon entangled and Fock states, whereas it is not 
a good quantum number in the coherent state. Finally, 
the third important aspect is entanglement that allows 
one to control the symmetry of the two-photon state. 
Symmetric and antisymmetric states that can be pre- 
pared in this way lead to coincidence rates R(t) that do 
not decrease with the pump bandwidth as fast as for the 
standard SPDC-generated entangled state. In addition, 
the antisymmetric state allows one to model fermionic 
behavior and change the constructive two-photon inter- 
ference into a destructive one. As a result, the peak of 
R{t) observed for the state with the bosonic symmetry 
(i.e., for the symmetric state), as well as for the states 
with no particular symmetry, at r = 0, turns into a deep 
for the state with the fermionic symmetry (i.e., for the 
antisymmetric state). 
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Appendix A: Photon number correlation versus 
coincidence rate 



In this Appendix, we establish a relation between the 
photon number correlation function Cy = (iUhj) and 
the probability P%(i, j) to detect a photon in each of the 
modes i and j. Let us first assume i ^ j. Denote the 
state of the field before detection of photons by \ip) and 
assume that {|?/>n)} is an orthonormal basis composed of 
all possible, orthogonal states in which the field can be 
found after the two photons are detected. The probabil- 
ity density to detect a photon from the mode i at a time 
t\ and a photon from the mode j afterwards, at a time 
t 2 >h, is then [14] 

P2(iJ;t 1 ,t 2 )=Y,\(i>n\& j (t2)a i (hM\ 2 

n 

= 5^<V|al(*i)a}(* a )|V'n>^n|a J -(t2)oi(ti)|^> 

n 

= (^|at(t 1 )at(t 2 )a J -(i 2 )a i (i 1 )|^) 

= (:n i (t 1 )h j (t 2 ):), (Al) 

where the operator cn{t) corresponds to the detection of 
a photon in the mode i at a time t [14] and we made use 
of the closure relation J2 n I 1 / 1 ™) (V^l = 1; tne c °l° ns : ' " " : 
denote normal ordering of operators. 

The probability that the two photons are detected at 
some times t\ < t 2 during a sampling time AT is ob- 
tained by integrating Eq. (Al) over times: 



AT/2 



AT/2 



r 2 



(M) 



J dh J dt2P 2 (hj;h,t 2 ). (A2) 



-AT/2 



Finally, the probability P2(i,j) of detecting the two pho- 
tons in arbitr ary order is equal to the sum of P* 1<t 2 (i, .j) 



given by Eq. ( A2 ) and P 2 £l > 2 given by Eq. ( A2 ) with 



the integration over t 2 running from — AT/2 to t\: 
P2(iJ) = P2 tl<t2 (iJ) + P2 1>t2 (iJ) 



AT/2 



AT/2 



J dti J dt 2 P 2 (i,j;ti,t 2 ) 

\T/2 ti 
YT/2 ti 

J dti J dt 2 P 2 (i,j;ti,t 2 ) 

\T/2 -AT/2 

\T/2 AT/2 

J dti J dt 2 P 2 (i,j;t 1 ,t 2 ) 



-AT/2 



-AT/2 

= (hihi 



(A3) 
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We therefore conclude that for i 7^ j, P2(i,j) and Cy are 
exactly equal. 

Assume now that i = j. We have 



where 



P 2 {i,i;t 1 ,t 2 ) = (: hiit^hiih) :) (A4) 



and 



AT/2 AT/2 

P 2 (M) = y d*i y dt 2 P2(i,i;t 1 ,t 2 ) 

-AT/2 ti 

= li--n*:) = l({n 2 i )-{fi i )). (A5) 

The difference with respect to the case i 7^ j comes from 
the fact that the two photons now belong to the same 
mode and hence are indistinguishable. The two cases 
t\ < t 2 and t\ > t 2 cannot be distinguished anymore and 
there is only one [instead of two in Eq. (A3 )] contribution 
to P 2 . Equations ( [A3] ) and |A5} lead to Eq. (flO} of the 
main text. 

The importance of the additional prefactor \ in Eq. 



( A5 ) as compared to Eq. ( A3 1 can be understood if we 



consider two identical one-photon wave packets incident 
at the same input channel of a symmetric beam splitter 
having two outgoing modes i, j = 1 (transmission) or 2 
(reflection) . Assume that the wave packets are well sepa- 
rated in time so that they interact with the beam splitter 
independently and each photon can be transmitted or re- 
flected with a probability Pi(l) = -Pi (2) = |. On the one 
hand, the calculation of joint probabilities readily yields 
P 2 {1,1) = £"2(2,2) = \ and P 2 (l,2) = §. On the other 
hand, we find (: n\ :) = (: ii 2 :) = (: hih 2 :) = |. We 
thus see that although P 2 (l,2) = (: fiih 2 :), an addi- 



tional factor \ 



is necessary to link P2(l, 1) and 



\\ :) 



^2(1,1) 



:). The difference between the cases 



i = j and i 7^ j comes from the fact that two different 
processes can lead to detecting one photon in the mode 
1 and the other — in the mode 2: either the first photon 
is transmitted and the second is reflected or vice versa. 
The (equal) probabilities of these two processes add up to 
give £2(1, 2). However, a unique process leads to finding 
both photons in the mode 1 (or 2): both photons should 
be transmitted (or reflected). 



Appendix B: Calculation of normally ordered 
photon number correlation functions 



In this Appendix, we derive Eqs. ( 11 )-( 13 ) of the main 
text. Consider first the case of t = 0, i.e. the case when 
there is no time delay between the two photon beams. 
We have 



AT/2 AT/2 

fiihj :) = y dti J dt 2 { 

-AT/2 -AT/2 



hi^hjih) :), (Bl) 



nSi)nj{t 2 ) :) = ($\a\ [h)a] [t 2 ) a,- (t 2 ) h (ti) |V) 



= \a {t 2 )a % {h)\^ 



(B2) 



Representing <ij(t) through its Fourier transform a i {io) 1 
hi(t) = (l/\/27r) f^°° du:ai{uj)e~ lut , we obtain 



aj(t 2 )ai(ti)\^) = 



2tt 



dudu aj(u))ai(tj ) 



x e 



-i(ut2+u'ti) I 



(B3) 



Using Eq. |l|, we now express dj(uj) and di(u>) through 
a (u) and a e (u)): 



dj{t 2 )di(ti)\ip) 



1 

2^ 



ckudw' 



x [t e j(u)a e (w) +t j(u>)a (u)] 
x [t cl {uj')d c (uj') + t oi (u)')a (u)')] 

x e -< wi2+a/tl )|^). (B4) 

Note that strictly speaking, we have to keep in the ex- 
pressions for dj(ui) and di(uj') in Eq. (B4) contributions 
from all the incoming modes others than the modes 'o' 
and 'e', even if no light is sent into these modes. Oth- 
erwise, the representations of the operators djiuS) and 
di(w') corresponding to outgoing modes, in terms of op- 
erators, corresponding to incoming modes, do not verify 
the bosonic commutation relations. We will, however, 
omit these additional terms from here on because they 
do not contribute to the normal-ordered correlation func- 
tions that we compute in this Appendix. 

For the rest of the calculation, it will be convenient to 
distinguish between the light in two-photon states (en- 
tangled and Fock states) and the light in the coherent 
state. 



1. Two-photon entangled and Fock states 

If the state |i/>) contains only two photons of which one 
has ordinary and the other one — extraordinary polar- 
ization, we can simplify Eq. (B4) by dropping the terms 



containing d {ui)d {ijj')\'ip) and d e {ui)d e {uj')\il)): 

POO 



aj{t 2 )ai(t{)\ip) = 



1 

2^ 



duduj 1 



x \t e j{uS)t oi {uj')d e {uj)d {u)') 
+ t e i(u)')t j(uj)aa(u)d e (Lj')} 
X e-'K^'* 1 )^), (B5) 

Using the bosonic commutation relations obeyed by the 
operators a c (oj) and d (u>'), we obtain 

d e (uj)a (uj')\ip) = // dujiduj 2 B(uji,uj 2 ) 
J J — 00 

x a c (u)a o (u')al(ui)al(u2)\0) 

= B(w , ,«)|0) J (B6) 
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and, similarly, 



a (w)a e (w')|^) =B(u,u')\0). 



(B7) 



We now assume that the sampling time AT during 
which the photons are counted is much longer than all 
other time scales of the problem. This allows us to take 



the limit AT — > oo in Eq. (Bl) and to use the integral 
relation dte^-^ 1 = 2M(u) -ft). We then obtain 



dw 1 dw 2 {\B(ui, lo 2 )\ 2 

' — oo 

x [|t OJ (w 2 )| 2 Mwi)| 2 + |i OJ (w 2 )| 2 |t C2 ( Wl )| 2 ] 

x [t ej {ui 2 )t* ei (u>i )t oi (u>i )t* OJ (ui 2 ) 

+ U^:>i)*o>i)4M}. (B8) 



We now average Eq. (B8) over disorder by using the 



Gaussian statistics of transmission coefficients t a i(uj), as 
discussed in Sec. |II A| This yields 



= 2T 



du\duj 2 B* (u)\, w 2 ) 



x B{u] 2 ,Wx)\C{ui - UJ 2 )\ 2 ] . 



(B9) 



The result for the Fock state is obtained by replacing 
B(uji,uj 2 ) by B(u>i)B(oj 2 ): 



n i n 3 VFock 



= 2T 



l + Si 



duj\dui 2 \B{ui)Y 



x |B( W2 )| 2 |C( Wl -w 2 )| 5 



(BIO) 



Equations (11) and (12 1 of the main text, correspond- 



ing to the situation in which the photon that has extraor- 
dinary polarization is delayed by a time r before entering 
into the disordered medium, are obtained by replacing 
t C j(uj) by i C j(w)e~ l " T in the above derivation. 



Coherent state 



Using the fact that the two-mode coherent state de- 
fined by Eq. ^ is an eigenstate of operators a 0iC (u;): 
a , e (a;)|V>coh) = a(uj)\ip coh } , we reduce Eq. ((B4j) to 



aj(t 2 )ai(ti)\il> C oh) = 



2tt 



diuduj' 



(Bll) 



t co {uj)t cl {uj')e-^+^ 

+ t Q] (u)t cl (J)e- lu ' T 
+ t j{u))t oi {u)')} 

x a(a;)a(w')e _4M2+ "'' tl) |V')- 

In contrast to the two-photon states considered in Sec. 
Bll the terms containing t e j(u))t e i(uj') and t j(ui)t i(ui') 



do not vanish. We now write down |aj(^2)ai(ii)|V'coh) 
integrate it over times t\ and t 2 in the limit AT — > 00, 
and average over disorder. This yields 



AT 2 d Wl dw 2 |a( Wl )| 2 |a( W2 )| 2 
J J — 00 

|2 „„„2/ 



x [l + <yc(wi -w 2 )| 2 cos 2 ([ Wl -wa]r/2)] . (B12) 



Finally, Eq. ( 13 ) of the main text is obtained after mak- 



ing use of the condition dui\a(ui)\ 2 — 1 that reflects 
the fact that each mode of our two-mode coherent state 
l^coh) contains one photon on average. 
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